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ABSTRACT 

We present line-of-sight stellar velocity distributions of elliptical galaxy NGC 821 obtained to ap- 
proximately 100" (over 2 effective radii) with long-slit spectroscopy from the Hobby-Eberly Telescope. 
Our measured stellar line-of-sight velocity distributions are larger than the planetary nebulae mea- 
surements at similar radii. We fit axisymmetric orbit-superposition models with a range of dark halo 
density profiles, including two-dimensional kinematics at smaller radii from SAURON data. Within 
our assumptions, the best-fitted model gives a total enclosed mass of 2.0 x 10 u Mq within 100", with 
an accuracy of 2%; this mass is equally divided between halo and stars. At lR e the best-fitted dark 
matter halo accounts for 13% of the total mass in the galaxy. This dark halo is inconsistent with 
previous claims of little to no dark matter halo in this galaxy from planetary nebula measurements. 
We find that a power- law dark halo with a slope 0.1 is the best-fitted model; both the no dark halo 
and NFW models are worse fits at a greater than 99% confidence level. NGC 821 does not appear to 
have the expected dark halo density profile. The internal moments of the stellar velocity distribution 
show that the model with no dark halo is radially anisotropic at small radii and tangentially isotropic 
at large radii, while the best-fitted halo models are slightly radially anisotropic at all radii. We test 
the potential effects of model smoothing and find that there are no effects on our results within the 
errors. Finally, we run models using the planetary nebula kinematics and assuming our best-fitted 
halos and find that the planetary nebulae require radial orbits throughout the galaxy. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 821) — galaxies: 
kinematics and dynamics — galaxies: halos — dark matter 



| 1. INTRODUCTION 

Cold dark matter is now accepted as an integral part of 
\ our universe, and recent observat ions have continued t o 
CN ■ provide support for its existence (|Komatsu et al.l 12008). 
^sO 1 Part of the picture of the universe is that galaxies are 
C*") ' surroun ded by massive dark matter halos in which the y 
^ ! formed (|White fc Reesl [l97l IBlumenthal eFaTl Il98l . 

■ Recently cosmological simulations have become detailed 
OO | enou gh to reach the level of individu al gal axy forma- 
O ' tion (|Naab et aLll2007t iGovernato et al.ll2007t) . and com- 
parisons with data can help further constr ain cosmo- 
logical theory (jOstriker fc SteinhardtJ 120031 ). Indeed, 
k> , spiral galaxy rotation curves are one of the strongest 
5_j ■ pieces of ob servational evidence fo r the existence of 
dark matter (Ivan Albada et al.l 119851: iPersic et ail 119961: 
ISofue fc Rubinll20'0lir It is also important to study the 
dark halo structure of elliptical galaxies because of their 
different formation and evolution. However it is more 
difficult to measure dark matter in elliptical galaxies be- 
cause of a lack of tracers at large radii where dark mat- 
ter is thought to dominate. The best way to measure 
the underlying gravitational potential is to use kinemat- 
ics from the stellar population, but this has been lim- 
ited due to the fai ntness of stellar ligh t in the outer 
regions of galaxies (|Gerhard et al.l 12001). Dark matter 
in elliptical galaxies has therefore been studied in other 
ways, such as via X-ray emission (|Loewcnstcin & White 

amydove@astro.as.utexas.edu, gebhardt@astro.as.utexas.edu 

1 Based on observations obtained with the Hobby-Eberly Tele- 
scope, which is a joint project of the University of Texas at 
Austin, the Pennsylvania State University, Stanford University, 
Ludwig-Maximilians-Universitat Miinchen, and Georg-August- 
Universitat Gottingen. 



119991: iMathews fc Brighentil I2003D. gravitat ional lens- 
ing (|Keetonl 120011: iMandelbaum et al.l 120081 ). and us- 
ing individu al stel l ar ve l ocities as in near by dwarf 
spheroidals (|Mateol 119981: IKlevna et all I2002D . In or- 
der to study a more representative sample of galax - 
ies, tracers such as globular clusters dZepf et al.l I200C 



Pierce et al.ll2006|) and planetary nebulae dMendez et al.l 
20011: IRomanowskv et all 120031: [Coccato et al.ll2009l) are 



used used to probe the outer parts of elliptical galax- 
ies, though it is difficult to get a significant sample size. 
Additional issues arise with these tracers, as discussed 
below, such as understanding their radial profile. With 
larger telescopes we are now able to measure stellar kine- 
matics from integrated light to larger radii, thus closing 
the gap between stars and the large-radii tracers. 

Meanwhile dynamical models of galaxies have also 
improved. Rather than previous spheric al models that 
use a nalytic distribution functions (DFs; iGerhard et al.l 
2001), orbit-based axisymmetric models are now avail- 
a ble. These fully gener al models, based on the technique 
of lSchwarzschildl (|1979f ) . provide detailed information on 
the orbital structure of the galaxy, including the DF and 
its projections, such as velocity anisotropy. Orbit-based 
models are now frequently applied to galaxies for studies 
of bo th dark halos and cent ral black ho les dRix et al 
19971 Ivan der Marel et all 119981: iCrettqn et al' 



Verolme et alJl2002tlGebhardt et al.l l2003l: iThomas et alJ 
2005: IGebhardt fc Thomasll2009D . 



19991 [Gebhardt et al.1 120001: iCappeilari et al.l 12002b 



The elliptical galaxy NGC 821 is an example in which 
the use of l arge-radii tracers h a s pro vided an intrigu- 
ing result. IRomanowskv et all (|2003t ) study the dark 
halo of NGC 821 using approximately 100 planetary neb- 
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ula velocities and find small line-of-sight velocity dis- 
persions tha t are c onsistent with little or no dark halo. 
iDekel et al.l ((2005) use disk galaxy merger simulations 
to show that large anisotropics can be created in the 
resulting elliptical galaxies, and that this anisotropy in 
combination with the different density profile of a young 
population could explain how the low dispersions from 
planetary nebulae measurements are also consistent with 
typical dark matter halos. Our study uses deep long- 
slit spectroscopy of NGC 821 from the 9.2-meter Hobby- 
Eberly Telescope to obtain stellar kinematics to greater 
than 2 effective radii in hopes of further constraining the 
d ark halo of this galaxy. 

iWeiimans et all (|2009fl model NGC 821 using data 
from SAURON, both at small radii (which we include 
in our analysis) and newer data at large radii. We find 
similar results both for the kinematics and for the dark 
halo properties. Comparison between the two studies is 

presented in their paper and withi n this paper. 

NGC 821 is classified as an E6? (Ide Vaucouleurs et al 



19911) . It has disky isophotes (|Lauerlll985l ; iBender et al 
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1988) and a power-law centr al surface brightness pro- 
file (|Ravindranath et al.ll200lD. The b lue absolute mag- 
nitude is -20.27 dTrager et all I2000D. We use a dis - 
tance of 23.44 Mp c taken from ICappellari et all (|2006f) . 
which adjusts the lTonrv et al.l (1200 ID values for the new 
Cepheids zero-point of lFreedman et al.l (120011). NGC 821 
is not detected in Ha (jMacchetto et al.1 119961 ) or OIII 
(jSarzi et al.ll2006h . Point source and diffuse X-ray emis- 
sion has been detected but there is no evidence for hot 
gas dPellegrini et all l2007allbT ). NGC 8 21 is considered 
a fas t rotator (ICappellari et al.l 120071 : lEmsellem et al.l 
12001 . iProctor et al.l (|2005t) find that NGC 821 has very 
strong age and metallicity gradients, from ^4 Gyr and 
3 times solar in the center to ~12 Gyr and less than ^ 
solar at lR e . They conclude that NGC 821 has experi- 
enced a recent (^1-4 Gyr ago) burst of star formation, 
most likely from in-situ gas and perhaps triggered by the 
accretion of a small satellite galaxy. This may be an in- 
dication that there are young planetary nebulae in this 
galaxy. 

f}2] describes the observations and data reduction; in 
$3] we describe the kinematic extraction; the dynamical 
models are described in 2J we present our results in $5] 
and give conclusions in [J6] 

2. OBSERVATIONS AND DATA REDUCTION 

Long-slit spe ctra were taken w ith the Low-Resolution 
Spectrograph (|Hill et al.l I1998D on the Hobby-Eberly 
Telescope. We use the g2 grism and 1" by 4' slit over 
the wavelength range 4300-7300A. This setup gives a 
resolving power of 1300 or a full-width half-maximum 
(FWHM) resolution of about 230 km s _1 . Measurements 
of night sky line widths show that we can measure dis- 
persions to about 110 km s _1 . The CCD frame (binned 
2x2) has a plate scale of 0.47"/pix spatially and 2A/pix 
spectrally. The gain is 1.832 e~ ADU -1 and readout 
noise is 5.10 e~. We use the Schott Glass blocking filter 
GG385, which has a half-power point of the transmission 
around 385 nm. 

NGC 821 was observed over eight nights in November 
2003 for a total exposure time of approximately 5.5 and 
2.3 hours on the major and minor axes respectively. Cad- 



mium and Neon calibration lamp exposures and white 
light illumination flat fields were taken each night. 

The data reduction uses standard techniques. First we 
overscan correct and trim the images. Then we apply a 
flat correction using a normalized flat frame, taken from 
averaged instrumental flats obtained each night of obser- 
vations. Next we rectify the images along the spatial axis 
using the calibration lamp lines as a reference. 

For sky subtraction, we use the region of the slit that 
is furthest from the galaxy center. Since we only have 
a 4' slit, there will be some galaxy light in the region 
where we select sky. However, the surface brightness pro- 
file extends out 350" so we can accurately calculate the 
amount of galaxy in our background region. For our last 
extracted spectrum (at 90"), the amount of galaxy light 
that we are including as background light is about 15% 
of the galaxy light for that last extraction. We have run 
simulations in order to determine whether this amount of 
contamination has an effect on the extracted kinematics. 
We take a high signal-to-noise galaxy spectrum and sub- 
tract off 15% of itself, and then extract the kinematics. 
Only for very high S/N does this amount have an effect 
and for the S/N for this dataset (as described below), we 
find no significant effect. 

3. KINEMATICS 

We extract the spectra in radial bins along the major 
and minor axes. Because the seeing is approximately 2" 
we set the central bins to 5 pixels (2.35"). The outer 
bins are sized to obtain sufficient signal for kinematic 
analysis. Along the minor axis the spectra from cither 
side of the galaxy were averaged at each radius before 
kinematic analysis. Along the major axis, the center of 
the galaxy was near the edge of the chip so only one 
side was extracted. Our farthest radial bin extends to 
99" on the majoraxis and 45" on the minor axis, corre- 
sponding to a V-band surface brightness of 21.8 and 23.5 
respectively. There is a faint object centered at 115" 
along the major axis that prevents a further radial bin 
from being used. The radial extent in effective radii de- 
pends on the value of R e used. The measured R e of 
NG C 821 varies throu ghout the literature: 50" (R C3) , 
45" dFaber et al]fl989h . 39" (ICappellari et al.ll2lfflfih . 36 " 
(jTrager et al.ll2000( ). and 16.7-18.3" (|Bender et al.l fl988). 
For the purpose of this discussion we adopt an R e of 45". 
Thus our data extend to approximately lR e along the 
minor axis and 2R e along the major axis. 

We do not flux calibrate the spectra, and thus we re- 
move the continuum in each spectru m. We fit the loca l 
continuum by finding the biweight (I Beers et al.l 119901 ) 
in windows as described in IPinknev et al.1 (|2003l ) . The 
wavelength solution comes from the Cd and Ne calibra- 
tion lamps. 

We obtain a nonparamctric linc-of-sight velocity dis- 
tribution (LOSVD) by deconvolving the galaxy spectrum 
with a set of stellar template spe ctra using the maximum 

^(goM). 

(pool . 



penalized likelihood technique of Gebhardt et al 



Tests of this technique are given in IPinknev et al 



There are 30 evenly-spaced velocity bins of 54 km s _1 
that represent the LOSVD. We vary the height in each 
bin and the weights of each template star to find the best 
match to the galaxy spectrum from each radial range. 
We use nine stellar te mplates with typ e s rang ing from G 
dwarf to M giant from lLeitherer et ail (|1996l ). convolved 
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Fig. 1. — Spectrum of the combined, weighted template stars 
(lower panel), data from the central bin along the minor axis 
(dashed line, upper panel), and the template spectrum convolved 
with the best-fitted LOSVD (solid line, upper panel). The region 
from 5163 A to 5228 A is excluded from the fit. 

to our spectral resolution. 

For our kinematic analysis we use the spectral range 
4800-5450A which matches the wavelength range of our 
template stars. This region includes the H/3 and Mgb 
lines, how ever we exclude the Mg6 region because it is 
enhanced (|Proctor et al.l 120051) and our template stars 
do not provide a proper fit. iBarth et al.l (|2002f) show 
that in pixel-space fitting routines the Mg6 line is sensi- 
tive to template mismatch and the details of the fitting 
procedure. If the Mg6 line is included in the fit, the mea- 
sured dispersions are falsely high by as much as 20% to 
account for the abundance discrepancy. An example fit 
is shown in Figure [1] 

The uncertainty of each velocity bin is obtained from 
Monte Carlo simulations. We convolve the best-fitted 
LOSVD and weighted stellar templates to obtain an ini- 
tial galaxy spectrum. We then generate 100 realizations 
of the galaxy spectrum by adding Gaussian noise using 
an estimate of the initial rms. The LOSVD is determined 
for each realization as described above. The distribution 
of values in each velocity bin of the LOSVD provides an 
estimate of the 68% confidence bands. The median value 
of the dispersion from the 100 realizations compared to 
the initial dispersion reveals any possible bias in the dis- 
persion measurement. 

Although we use the full nonparametric velocity profile 
in the dynamic modeling it is useful to compare moments 
of the distribution. In Figure [2] we plot, from top to bot- 
tom, the second moment as measured by \f[V 2 + cr 2 ), 
the first four Gauss-Hermite moments (mean velocity V, 
velocity dispersion a, asymmetric deviations from Gaus- 
sian (similar to skewness) /13, and symmetric deviations 
from Gaussian (similar to kurtosis) /14). The kinematic 
data are given in Table [T] and Tab led For comparison , 
in F igure [2] we also p l ot da ta from lPinknev et all ()2003[ ) 
and lEmsellem et al.l (|2004| ) extracted in a 1 arcsec slit 
along the major and minor axes. The second moment of 
the line-of-sight velocity, (V 2 + C7 2 ) 1 / 2 , is slightly smaller 
than the other samples throughout the overlapping re- 
gion. This may be caused by a slit misalignment (since 
V will be higher on the major axis) or template fitting 
difference. Since, however, our models are dominated 
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Fig. 2. — Gauss-Hermite moments of the LOSVDs and the rms 
line-of-sight velocity: mean velocity V, velocity dispersion tr, asym- 
metric deviations from Gaussian (skewness) h^, and symmetric de- 
viations from Gaussian (kurtosis) /14 along the major axis (left 
panel) and mino r axis (right panel) fo r our data (black filled cir- 
cles), SAURON (Emscllcm ct al. 2004) (blue open triangles), and 
Pinkncv ct al. (2003) (red open squares). 

by SAURON data in the center, this difference is not a 
major issue. We further run the dynamical models with 
using SAURON data alone and HET data alone, and find 
un-biased results from when using the combined dataset. 

Our hi values are more negative than the data from 
the literature. This discrepancy could be attributed to 
template mismatch (either in the published analysis or 
in ours), however we use a wide range of template stars 
and do not get a different result when more template 
stars are made available for the fit. It could also be 
that relying on Gauss-Hermite parameterization causes 
some differences since there are known correlations, es- 
pecially with higher o rder moments (seefMagorrian 2006; 
iHoughton et all 12006( 1. Since we fit the LOSVD directly 
in the dynamical models, a better comparison would be 
with those profiles, as opposed to their moments. The 
dark halo mass, however, is determined mainly by the 
radial profile of the second moment, and /14 determines 
mainly the anisotropics. There is certainly some degener- 
acy between the two parameters, but we find no reason 
to believe that our hi values are incorrect. There are 
also kinematic points that have differences which are in- 
consistent with their reported uncertainties (for example, 
some of the minor axis points), and the uncertainties may 
be underestimated for those points. We run halo mod- 
els without the most discrepant points and still find the 
s ame halo result s as wh en they are included. 

iProctor et all ([2009) measure the kinematics of 
NGC 821 using the background galaxy light from Keck 
DEIMOS multi-object spectroscopy of globular clusters. 
Their results show good agreement with ou r V and a pro- 
files, a s seen in their Figure 17. Previously. IProctor et al.l 
( 2005) determined the kinematics of NGC 821 using the 
Gemini Multi-Object Spectrograph (GMOS). As shown 
in Figure 2 of that paper, our a is in agreement with 
theirs, except that they find a lower dispersion at about 
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TABLE 1 

Major Axis Kinematics 



r 


V 




cr 


Co- 


h 3 


"-3 


tl4 




arcsec 


km s~ 




km s~ 












0.00 


-1.05 


1.65 


208.24 


4.43 


0.049 


0.040 


-0.076 


0.012 


1.41 


25.35 


1.40 


191.71 


2.96 


0.005 


0.025 


-0.056 


0.010 


3.76 


48.53 


1.55 


185.06 


4.16 


0.002 


0.038 


-0.061 


0.008 


6.11 


60.13 


1.89 


186.95 


4.73 


-0.007 


0.030 


-0.063 


0.008 


8.46 


65.68 


1.72 


177.81 


3.93 


-0.010 


0.016 


-0.035 


0.009 


10.81 


60.16 


2.79 


179.86 


4.20 


-0.023 


0.018 


-0.042 


0.008 


13.16 


85.90 


3.80 


172.81 


4.40 


0.008 


0.016 


-0.047 


0.010 


15.51 


73.79 


3.76 


175.29 


5.35 


-0.005 


0.018 


-0.051 


0.010 


18.09 


61.67 


3.37 


184.25 


6.40 


-0.002 


0.025 


-0.035 


0.012 


21.62 


77.70 


3.91 


172.87 


5.93 


0.004 


0.020 


-0.058 


0.011 


26.08 


56.48 


3.27 


177.18 


4.56 


-0.016 


0.021 


-0.044 


0.009 


31.96 


55.20 


3.46 


171.32 


5.12 


-0.022 


0.021 


-0.034 


0.013 


39.01 


42.58 


4.29 


168.27 


4.21 


-0.039 


0.019 


-0.039 


0.007 


47.24 


66.29 


3.89 


160.18 


5.64 


0.068 


0.015 


-0.040 


0.009 


58.99 


42.41 


6.68 


176.27 


6.69 


-0.012 


0.017 


-0.041 


0.010 


74.26 


33.40 


8.50 


173.74 


6.79 


0.100 


0.022 


-0.003 


0.018 


90.47 


79.97 


6.16 


170.17 


7.48 


-0.019 


0.026 


-0.019 


0.017 



TABLE 2 
Minor Axis Kinematics 



r 


V 




a 




h 3 


e h 3 


fl4 


£h4 


arcsec 


km s~ 1 




km s — 1 












0.00 


-6.51 


1.29 


211.91 


5.50 


0.069 


0.053 


-0.062 


0.016 


1.41 


-4.68 


0.86 


201.49 


4.47 


0.037 


0.022 


-0.066 


0.009 


3.76 


-7.94 


1.04 


188.69 


3.30 


0.000 


0.015 


-0.053 


0.009 


6.11 


-6.27 


2.88 


193.16 


4.54 


-0.012 


0.024 


-0.053 


0.011 


8.46 


1.79 


2.63 


198.46 


4.27 


0.005 


0.027 


-0.049 


0.011 


10.81 


-6.24 


2.72 


185.65 


3.66 


-0.015 


0.020 


-0.046 


0.011 


13.16 


-1.07 


4.47 


195.26 


6.04 


0.083 


0.029 


0.025 


0.024 


15.51 


14.06 


5.22 


209.36 


5.76 


-0.033 


0.028 


0.010 


0.020 


17.86 


6.49 


6.04 


180.82 


7.07 


-0.043 


0.033 


-0.059 


0.016 


20.21 


-31.59 


6.41 


182.22 


12.62 


0.055 


0.045 


0.021 


0.024 


22.56 


30.06 


12.08 


219.63 


11.81 


-0.031 


0.037 


-0.030 


0.026 


26.08 


-40.85 


10.30 


177.36 


11.96 


-0.044 


0.035 


-0.021 


0.022 


30.78 


5.47 


10.77 


194.77 


12.27 


-0.036 


0.037 


-0.034 


0.024 


35.49 


-23.16 


14.05 


171.31 


11.39 


0.100 


0.041 


0.016 


0.028 


41.36 


-76.04 


21.58 


233.45 


21.95 


0.020 


0.050 


-0.004 


0.032 



30". However their data are of substantially lower signal 
to noise and they are also unable to determine higher 
o rder moments. 

iWeiimans et al.l (|2009f ) provide a new analysis of the 
SAURON data and also include additional data at large 
radii. Their furthest radial point is at 110" (which they 
refer to as 4 Re), whereas our last point is at 90" (which 
we refer to as 2 Re). The comparison between the two 
kinem atic sets is shown in Figure 6 from IWcij man s et al.l 
(2009). There is generally excellent agreement between 
the two sets of kinematics. Furthermore, their re-analysis 
of the SAURON central pointing shows /14 values now 
more consistent with our numbers. The higher-order 
moments of the LOSVD are difficult to measure, and 
it is important to consider systematic diffe r ence in the 
analysis. The spectra from IWeiimans et~a l. (2009) have 
lower signal-to-noise than our spectra, which could add 
to systematic difference. For this reason, our dynamical 
modeling does not include their kinematics, although we 
suspect there will be little difference in the overall results. 
Thus, we can compare the constraints on the dark matter 
parameters, which would include systematic differences 
in the kinematic samples used. 



Our results are consistent w ithlWeiimans et al.l (12009) . 
an d not consisten t with either IRomanowskv et all ()2003T ) 
or lCoccato et al.l (|2009l ). The Weijmans et al. data are 
consistent with that of Coccato et al. The reason is sim- 
ply that the Weijmans et al. dispersions have large un- 
certainties, and that the Weijmans et al. dispersions are 
between our values and that of Coccato et al. The differ- 
ence is not due to comparing dispersions measured from 
major axis radii compared to circular radii (as in Coc- 
cato et al.), since our dispersion profile is nearly flat. We 
compare with the Romanowsky et al. values and not to 
the values reported in Coccato. The differences between 
Coccato et al. and Romanowsky et al. are small enough 
to not impact our analysis or conclusions. The Coccato 
et al. data have slightly smaller uncertainties than Ro- 
manowsky et al., so the statistical differences between 
our dispersions and theirs is slightly larger. 

Figure [3] shows the rms line-of-sight velocity (V 2 + 
a 2 ) 1 / 2 compared to the planetary nebula results of 
IRomanowskv et ail (|2003f) . Our largest radii data show 
higher rms line-of-sight velocities than the planetary neb- 
ulae, at about 3cx for their two largest radii points. Thus, 
there appears to be a significant difference in the kine- 
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TABLE 3 
V-band Luminosity 
Density 

radius (") L Q /pc 3 



2.300E-02 
2.533E-02 
2.790E-02 
3.073E-02 
3.385E-02 
3.728E-02 
4.107E-02 
4.523E-02 
4.982E-02 
5.487E-02 



8.275E+03 
7.205E+03 
6.283E+03 
5.485E+03 
4.793E+03 
4.192E+03 
3.668E+03 
3.212E+03 
2.813E+03 
2.464E+03 



Fig. 3. — The rms line-of-sight velocity (v 2 + tr 2 ) 1 / 2 in km s — 1 
as a function of radius from our major axis data (filled squares) 
and minor axis data (open squar es), compared to the data from 
planetary nebulae measurements (Romanowsky ct al. 2003). Wc 
use v and a as measure from a Gauss-Hermite fit; since we do not 
correct for the higher order moments, these values approximate the 
actual second moment. The arrow indicates the adopted R e of the 
galaxy. 

matics between the two samples. 

4. DYNAMICAL MODELS 

We use axisymm etric orbit superpositio n models based 
on the method of Schwarzschild (1979). The surface 
brightness profile is converted to a luminosity density 
profile using an assumed inclination. We assume an edge- 
on inclination for this analysis, which is reasonable given 
NGC 821's large ellipticity of 0.40 (|Cappellari et al.l 
l2007f h This luminosity density is converted to a mass 
density using a mass-to- light ratio (M/Ly) that is con- 
stant over the galaxy. A spherically symmetric dark halo 
density profile is added to the stellar density and this to- 
tal mass density gives the galaxy's gravitational poten- 
tial. Next individual stellar orbits are sampled in energy 
(E), angular momentum (L z ) and the third integral (I3) 
and these orbits are integrated in the specified poten- 
tials. The galaxy is divided spatially into cells both in 
real space and in projection, and the amount of time that 
an orbit spends in a cell represents the mass contributed 
by that orbit. The orbits are combined with nonnegative 
weights to find the best-fitted superposition to match 
the data LOSVDs from both HET and SAURON and 
the light profile. The model incorporates seeing by con- 
volving the light distribution for every orbit with the 
appropriate PSF before comparing with the data. This 
process is repeated for different dark halo density profiles 
and Mj Ly values to find the halo potential that best fits 
the data, as determined by \ 2 (describe d in fl5l). 

The SAURON data are described in lEmsellem et al.1 
( 2004). We reconstruct a full LOSVD from their reported 
moments and r ebin t hem to match our model bins. The 
iPinknev et alj (|2003|) data shown in Figure [2] is not used 
in the models. 

To reduce computational time, an orbit library is cal- 
culated for a given input dark halo plus stars with a 
mass-to-light ratio of one. The velocities are then scaled 
accordingly given the mass-to-light ratio before match- 
ing the data. The numbers reported are the actual den- 
sity parameters, including this M/L factor, which gives 



the somewhat irregular parameter space grids (as seen in 
Figure [5]). 

We u se th e orbital weight fitting of iGebhardt et al.1 
(2000. 120031) with the orbit library sampling of 
IThomas et all (12004120051). Our models differ from oth- 



ers (e.g. Cretton et al.l 119991) in that we use maximum 
entropy (jRichstone fc Tremaind n-988) an d we utilize the 
full LOSVD , rather than its moments. IThomas et al.1 
(|2004l I2005D show the ability of our orbit libraries to re- 
covery dark halo profile from mock elliptical galaxy data. 
Therefore these models should accurately measure the 
properties of NGC 821 given the caveats that we assume 
an axisymmetric galaxy and spherical halo. 

The orbits are computed in 4 angular bins and 15 ra- 
dial bins from 0.3" to 300". These bins are similar in 
size to the HET data extraction bins in the radial direc- 
tion, though in the angular direction they span about 20 
degrees. We specify the galaxy potential and the forces 
on a grid that is four times finer. Our libraries have 
approximately 10000 total orbits. 

To calculate our galaxy potential we use a composite 
surface brightness profile. Within 0.3" we use the pro- 
file fro m jLauer et al.1 (|2005l ) as compiled in lPinknev et al.1 
(2003) based on HST WFPC2 images in F555W. Outside 
of 0.3" we use a composite profile from HST PC F555W 
and the McDonald Observatory 0.8-m telescope in V (D. 
Fisher, private communication). The surface brightness 
deprojection is based on a nonparame tric estimate of 
the d ensity using smoothing splines (see fGcbhard t et all 
1996). The luminosity density is given in Table [3] 

4.1. NFW Halo 

We use the NFW (jNavarro et al.lll996| ) dark halo den- 
sity profile, given as 

Pcrit &c 



p(r) 



(r/r s )(l + r/r s ) 2 



(1) 



where r s is the scale radius of the halo, p cr n = 3H 2 /8irG 
is the critical density, and 5 C is the characteristic over- 
density. We use H = 70 km s _1 Mpc -1 . Throughout this 
paper we refer to p cr it$c &s the scale density. We use two 
independent parameters to define the NFW halo in our 
moedls: the scale density and the scale radius. The scale 
density can also be written in terms of a concentration 
parameter c by 



(2) 



3 ln(l + c) -c/(l + c)' 



() 
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The virial overdensity A V i r varies with redshift and 
cosmological model and we use a value of A v i r = 
101. Although we vary both the concentration (density) 
and sc ale radius, there is a known correlation between 
them jN avarro et all 11996(1 . This relation as given m 
IBullock et al.l (|2001l) is 



M„ 



1.5 x IO^H^Mq 
and can be written in the form 



-0.13 



(3) 



-1/0.13 



47T 



A mr — Pcrit c 3 ) (1.5 x 10 13 h- x M Q ). 

(4) 



4.2. Power-Law Halo 

The best-fitted NFW halo profiles have a break radius 
beyond the extent of our modeling and therefore look 
like a power-law over the extent of our models (see Ej5.ll 
below) . We therefore tried a simple power-law profile as 
well. We used power-law density profiles of the form 



P(r) = Po(—)~ 



(5) 



where n is the power-law slope, p is the characteris- 
tic density, and r Q is the characteristic radius such that 
p(r = t ) = p - We use r D = 0.3" = 34 pc because it is 
the inner-most radial point calculated in the models. 

5. RESULTS 

The best-fitted model is determined by comparing the 
X 2 between the model and data LOSVDs, with the uncer- 
tainty of the data determined from the 68% confidence 
band. Example LOSVDs are shown in Figure [4] for sev- 
eral radial bins. The measure of the reduced x 2 is not 
straight-forward since because determining the number 
of degrees of freedom is uncertain. The number of in- 
dependent observables is roughly the number of radial 
data bins times the number of LOSVD bins at each ra- 
dius (69 x 13 = 897 in this case), however the LOSVD 
bins are correlated and thus the effective number of data 
points is less than this value. The best-fitted model has 
a x 2 of around 2200, and with 897 data points, this pro- 
vides a large reduced x 2 ■ Typical values of the reduced 
X 2 for the orbit-based models are around 0.5 (see Geb- 
hardt et al. 2003), so the value reported here is not 
typical. The main driver for the large x 2 is the minor 
axis data — removing this data gives a reduced x 2 be- 
low one. Furthermore, the results on the parameters 
do not change significantly. Regardless, the change in 
X 2 between different models remains a valid statistic to 
determine confidence levels of the fits. For example, a 
change in x 2 of 2.3 corresponds to the 68.3% confidence 
level because we marginalize over M/L and thus have 
two parameters describing the halo. 

Because of computational limits we first calculate mod- 
els using a coarse grid of mass-to-light ratio. The x 2 
values are then fit with the IDL quadratic interpolation 
routine, and those models with the lowest minimum x 2 
are modeled with a finer mass-to-light ratio interval. 

5.1. NFW Halo 
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Fig. 4. — Match of data and no dark halo model LOSVDs for the 
central four radial bins along the major axis (r = 0.00, 1.41, 3.76, 
6.11 arcsec). The open circles are the data values with error bars 
and the closed circles are the model values. The area is normalized 
to the total light in that bin. 
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Fig. 5. — Scale radius (r 3 ) and scale density (p C rit$c) X 2 grid for 
NFW halo density profiles. Each point represents a model, and the 
size of the point reflects the value of A\ 2 for the best-fitted M/Ly 
value. Models with A\ 2 less than 6<r from the minimum value are 
plotted with open circles. The ringed point indicates the model 
with the lowest value of \ 2 - The d ashed line shows the expected 
NFW parameter relation (see ^4~f}. 

We use models with scale radius from 1 to 2000 kpc 
and scale density from 0.05 to 3.0 x 10~ 5 M Q pc _3 , cor- 
responding to a range in c of approximately 0.75 to 23, 
and MjLy from 1.0 to 9.0. Figures [5] and [6] show the 
resulting x 2 a s a function of halo scale radius and scale 
density. The points represent actual modeled values, and 
the M/Ly that gives the lowest x 2 is used at each point. 
The dashed line in Figure indicates the expected cor- 
relation of concentration and scale radius as described 
in t j4.ll This relati on has a scatter of Alogr s = 0.36 
(Bullock et al. 2001). Our data show a degeneracy be- 
tween scale radius and scale density that is similar to, 
though slightly tilted from, the correlation. 

X 2 is a function of 3 variables: stellar M/L, dark halo 
scale radius and dark halo normalization. Due to com- 
puter resources, we do not provide a uniformally-sampled 
grid of the 3 variables for the x 2 ■ Because of this, it is 
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Fig. 6. — x 2 as a function of scale radius and scale density 
(PcritSc) for NFW halo density profiles with best-fitted M/Ly. 
The points represent actual modeled values. The dashed lines refer 
to Ax 2 = 2.3, 4.61, 6.17, 9.21, 11.8, 18.4, corresponding to 2 degree 
of freedom confidence levels of 63.8%, 90%, 95.4%, 99%, 99.73%, 
and 99.99%. 
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difficult to produce reliable contours for any 2 of the pa- 
rameters. Figure [5] thus shows only the location of the 
points (with size related to x 2 )- We do not estimate un- 
certainties from the contours directly, but instead rely 
on plotting x 2 versus each of the parameters, including 
all values for the other two parameters. Figure [5] shows 
X 2 versus scale radius and density. Uncertainties come 
from the envelope of these one-dimensional plots. Since 
we have explored neither a regular grid nor a full set of 
variables (e.g., black hole mass, inclination, change in the 
stellar M/L with radius), the uncertainties should only 
be used in a comparative sense with the models that we 
have tried. A full exploration of the uncertainties will 
come as computer resources improve. 

We find that the best-fitted NFW dark halo density 
profile has scale radius 700^3qq kpc and scale density 
1.28l°;g x 10~ 4 M©pc" 3 , corresponding to a c of 2.45. 
The no halo model is ruled out with a change in x 2 of 
356 (greater than 99% confidence level) from the best- 
fitted NFW halo. Table E] shows the x 2 values and halo 
parameters of the best-fitted halo model and model with 
no dark halo. We do not attach significance to this den- 
sity, radius, and concentration. They are clearly outside 
the expected range for a galaxy and merely indicate that 
the NFW profile is not reasonable. The halo needs more 
mass at large radii to fit the data. The best-fitted NFW 
halo density profile is shown in FigureO The scale radius 
is well beyond the radial extent of our modeling, and is 
indicative of the need for a near power-law profile over 
the extent of our models. 

Since dynamical modeling directly measures mass (as 
opposed to dark halo parameters) , the enclosed mass pro- 
vides a more robust estimate and is likely not subject 
to the specific parameterization of the dark halo. Fig- 
ure |S] shows the mass enclosed within the extent of our 
kinematic data as a function of \ 2 . The best-fitted to- 
tal enclosed mass is 1.78 ± 0.15 x 10 11 M Q , divided into 



Fig. 7. — Density (top) and circular velocity (bottom) as a func- 
tion of radius for the best-fitted NFW (dashed lines) and power-law 
(solid lines) dark halos. In each case the bottom line is the dark 
halo alone and the top line is the total mass (halo plus stars) The 
data point shows the radius of the most extended bin of our kine- 
matic data. 
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total mass (M Q ) mass in stars (M Q ) mass in halo (M e ) 

Fig. 8. — Enclosed mass within the radial extent of our kine- 
matic data, 100", as a function of \ 2 f° r both the NFW halo 
density profiles (top) and power-law halo density profiles (bot- 
tom) with best- fitted M/Ly. The dashed lines refer to Ax 2 = 
2.3, 4.61, 6.17, 9.21, 11.8, 18.4, corresponding to 2 degree of freedom 
confidence levels of 63.8%, 90%, 95.4%, 99%, 99.73%, and 99.99%. 

1.03 ±0.03 x 10 U M Q in stars and 0.75 ±0.15 x 10 n M o 
in dark matter. At lR e the ratio of dark matter to to- 
tal matter is 0.19. The best-fitted NFW halo circular 
velocity profile is shown in Figure [7] 

Figures [9] and [TU] show the internal moments a r , ag , 
and 00 and ratio of radial to tangential dispersion along 
the major and minor axes for the model with no dark 
halo and the best-fitted NFW halo model. The model 
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Fig. 9. — Internal moments ay, <rg, and cr^, (top) and the ratio of 
radial to tangential dispersion (bottom) along the major axis (left) 
and minor axis (right) for the model with no dark halo. Note that 
a$ includes both random and ordered motions, which are shown 
(dot-dashed line) and are small. The vertical line shows the limit 
of our kinematic data; results beyond this radius are not reliable. 




Fig. 10. — Internal moments ay, erg, and (top) and the ratio of 
radial to tangential dispersion (bottom) along the major axis (left) 
and minor axis (right) for the model with the best-fitted NFW 
halo. Note that includes both random and ordered motions, 
which are shown (dot-dashed line) and are small. The vertical line 
shows the limit of our kinematic data; results beyond this radius 
are not reliable. 

without a dark halo shows radial anisotropy at small radii 
and tangential anisotropy at large radii along the major 
axis. Tangential anisotropy at large radii in a model with 
no dark halo could be an indication of the need for a dark 
halo because the observations largely constrain only 
(for an edge-on configuration), so both oy and o~q may 
be artifically decreased to create a smaller total a that 
can be fit without a dark halo. Along the minor axis, the 
contribution in the 9 and <j> directions are roughly equal, 
as is expected for an axisymmetric model. Overall the 
minor axis shows tangential anisotropy over the entire 
range of our data. The NFW model is more isotropic 
than the model with no halo. 

Although the models fit the full nonparametric veloc- 
ity profile of both the HET data and SAURON data, in 
Figure HT1 we plot the first four Gauss- Hermite moments 
for our HET data and the best-fitted halo models. The 
models differ most at intermediate to large radii, and do 




Fig. 11. — Gauss-Hermite moments (mean velocity V, veloc- 
ity dispersion a, asymmetric deviations from Gaussian (skewness) 
hs, and symmetric deviations from Gaussian (kurtosis) /14) of the 
LOSVDs for our HET data and the best-fitted halo models along 
the major axis (left panel) and minor axis (right panel). The HET 
data are shown with open circles, SAURON data along the axes 
with open triangles, no dark halo model with solid lines, best-fitted 
NFW halo with dashed lines, and best-fitted power-law halo with 
dotted lines. The model fits the full LOSVD of the HET data and 
all of the SAURON data. 

not appear to be driven by any one single parameter or 
radius in particular. 

5.2. Power- Law Halo 

We run models with a range of slope n from 0.0 to 1.1, 
density p Q from 0.0015 to 26 M /pc 3 , and M/Ly from 
3.5 to 8.0. The resulting x 2 grid is shown in Figure IT21 
and as a function of n and p a in Figure H31 The best- 
fitted halo model has a slope O.llgjg and a characteristic 
density p Q = 0. 025±° ;So|M©pc -3 . This halo is a better 
fit to the data than the best NFW halo, with a Ax 2 = 
95 (see Table 2]) . This power-law slope is significantly 
more shallow than the 1.0 slope of an NFW profile. A 
comparison of the best-fitted halo density and circular 
velocity profiles is shown in Figure [7j 

The best-fitted total enclosed mass is 2.01 ± 0.15 x 
10 n M o , divided into 1.04 ±0.02 x 10 n M Q in stars and 
0.97 ± 0.15 x lO n M in dark matter (see Figure HJ). At 
lR e the ratio of dark matter to total matter is 0.13. The 
internal moments a r , o~g, and cr^ and ratio of radial to 
tangential dispersion along the major axis are shown in 
Figure [14] and are roughly consistent with those of the 
best-fitted NFW halo. 

Figure [11] shows the first four Gauss-Hermite moments 
for our HET data and the best-fitted halo models. Note 
that the models fit the full nonparametric velocity profile 
of both the HET data and SAURON data. 

5.3. Model Tests 

In order to learn which aspect of the data is driving 
the results we have performed several tests. Using an 
abbreviated grid of about one third of the halos in the 
full run and a coarse spacing in M/Ly we have re-run 
the models with various subsets of the data. First, to 
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TABLE 4 
Best-Fitted Halo Model Results 



halo 
(1) 


x 2 

(2) 


M/Ly 
(M/L) e 
(3) 


r s 
kpc 
(4) 


p 

M Q /pc 3 
(5) 


c 

(6) 


1O 17 M 
(7) 


n 
(8) 


Mtot 
lO n M 

(9) 


Mstars 

1O 11 M 
(10) 


Mhalo 

lO n M 


none 
NFW 
power-law 


2527.30 
2171.70 
2077.05 


7.25 ±0.05 
6.19 ±0.09 
6.25 ±0.07 


70Q +500 


1.28 + g'| x 10~ 4 
025+ 025 


o 45+0.65 


4.76±f 


1+ 01 
u,1 -0.08 


1.20 ±0.01 
1.78 ±0.15 
2.01 ±0.15 


1.20 ±0.01 

1.03 ±0.03 

1.04 ±0.02 


0.00 
0.75 ± 0.15 
0.97 ±0.15 



Note. — (1) Dark halo density profile. (2) x 2 of best-fitted model. (3) Stellar M/Ly of best-fitted model. (4) Scale radius of best-fitted 
model. (5) Scale density p cr it5 c for NFW, characteristic density p a for power-law. (6) NFW concentration parameter determined from scale 
density. (7) Virial mass determined from NFW concentration parameter. (8) Power-law index. (9) Total mass within 100". (10) Mass of 
stars within 100". (11) Mass of dark halo within 100". 
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log Density (M gun pc ) 

Fig. 12. — Power-law slope n and density p a x 2 grid for power- 
law halo density profiles. Each point represents a model, and the 
size of the point reflects the value of Ax for the best-fitted M/Ly 
value. Models with Ax 2 less than 6<r from the minimum value are 
plotted with open circles. The ringed point indicates the model 
with the lowest value of x 2 - 

address any concerns over the scattered minor axis data 
we have run the test models using only the HET major 
axis data and SAURON data. The results are the same 
as the full data set; there is a clear need for a dark halo, 
and the best-fitted NFW halo is not as good a fit as the 
power-law halo. Second, we removed the two points on 
the major axis with extreme /13 values (at about 47" and 
74") since Figure [TT] may lead one to believe they are 
driving the fits. Again the results are the same as with 
the full data set. And third, we do a test run using only 
data below 0.5i? e . In this case there is essentially no 
difference in x 2 between the three best-fitted models (no 
dark halo, NFW halo, and power-law halo), and the best- 
fitted halos are quite different than those from the full 
data results. These tests indicate that it is the large radii 
data as a whole that is driving the model fits. To further 
demonstrate this, Figure [15] shows the A% 2 between the 
model with no halo and the best-fitted halo in each bin. 
The bins at large radii show the greatest change in x 2 , 
again indicating that the large radii data are the major 
factor in the fits. 

5.4. Comparisons to Other Studies 

iGebhardt et al.l (|2003l) model the central region of 
NGC 821 and find that it is radially anisotropic within 
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Fig. 13. — x as a function of power-law index n and scale density 
p for power-law halo density profiles with best-fitted M/Ly. The 
points represent actual modeled values. The dashed lines refer to 
Ax 2 = 2.3,4.61,6.17,9.21, 11.8, 18.4, corresponding to 2 degree of 
freedom confidence levels of 63.8%, 90%, 95.4%, 99%, 99.73%, and 
99.99%. 




Fig. 14. — As in Figure [9] for the best-fitted power-law model. 

a few arcseconds and isotropic to slightly tangentially 
anisotropic at larger radii. Given the difference in 
spatial resolution this roughly agrees with our result. 
iCappellari et al.l (|2007[ ) find that within about 20" N821 
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Fig. 15. — Difference in x 2 between the LOSVDs of the model 
with no dark halo and the best-fitted power law halo model aver- 
aged in each spatial bin. Green indicates that the no-halo model 
has a larger \ 2 than the power-law model and therefore the power- 
law is a better fit, while orange indicates that the power-law model 
has a larger \ 2 than the no-halo model and therefore the no halo 
model is a better fit. 

is radially anisotropic overall, along the major and minor 
axes and between. They find that the velocity ellipsoids 
are circular in the center and become more radial with 
increasing radius, in conflict with our results, but they 
do not include a dark halo which could change their re- 
sults. I Thomas et a l. (2007), using similar modeling as 
we use, find that early- type galaxies in the Coma cluster 
are radial compared to the 9 direction over all radii along 
the major axis, agreeing with our result, and vary from 
galaxy to galaxy in the <j> co mponent. The merger simu- 
lations of lDekel et a l. (2005) also find a radial anisotropy. 
Their spherically averaged j3 of about 0.4 corresponds to 
a (Jradiai / (Ttangentiai of 1.3, which is larger than our re- 
sults along the major and minor axes. However their 
simulations show declining projected dispersion profiles, 
which our data does not, that could account for the dif- 
ference. 

By modeling only the central part of NGC 821 
iGebhardt etaLl ([20031 ) find M/Ly = 7.6 (without in- 
cluding foreground extinction), which is consistent with 
our no-halo value over the whole galaxy of M/Ly = 7.25. 
Corre cting for NGC 821 's large reddening of Ay = 0.364 
mag (|Schlegel et al J 119981 NED extragalactic database) 
we find our best- fitted M/Lv nn hain = 5.18, M/Ly„,f w = 
4.43, and M/L ViPOW = 4.47. iCappellari etall (2006) 
find M/L jeans = 3.54, M / L schwarzschad = 3.08, and 
M/L steiinmcm = 2.60 in the / band. Using (V— J) = 1.35 
mag (jLauer et all 120051 ). deredden ed to (V - J) n = 1.20 
mag, and (V — /)© = 0.682 mag (jRamirez fc Melended 
2005) our V band mass-to-light ratios are converted to 
M/L Ijnohalo = 3.21, M/L I:Tlfw = 2.75, and M/L I>pow = 
2.77. Our mass-to-light ratios are slightly higher, though 
roughly consistent, with their mass-to-light ratios found 
using Schwarzschild modeling and stellar populations. 

We find that the enclosed mass of NGC 821 within 
~ 2R e is roughly 2 x 10 1 1 Mq, equally divided between 
stars and dark matter. At lR e the ratio of dark matter 
halo mass to total mass is 0.19 for the best- fitted NFW 
halo profile and 0.13 for the best- fitted power-law halo 
profile. This matches other studies that find that the 



dark matter is 10 — 40% of the total matter at lR e and 
that dark matter begins to dominate at 2 — AR e (e.g. 
Sariia et al.ll2000t IGerhard et alJl200lt lM amon fc Lokad 
20051 ). The simulations of iDekel etahl (|2005| ) also show 
that dark matter and stellar matter are equal at 3i? e , 
and at 1R P _ h a ve a mass fraction of 40% dark matter. 
I Thomas et al.l ()2007f ) perform similar dynamical model- 
ing on 17 galaxies in the Coma cluster. Using values 
taken by eye from their Figure 5 we find that their av- 
erage dark matter fraction at lR e is 0.19, though their 
galaxies show a wide range of dark matter fractions, from 
about 0.1 to 0.5 at li? e . We therefore find that the dark 
matter fraction at lR e is similar for N821, a field ellip- 
tical galaxy, and a selection of Coma cluster early-type 
galaxies, perhaps contrary to hypotheses that environ- 
ment plays a role i n the dark matter fraction. 

iWeiimans et al.1 ()2009f ) provide a dynamical analysis 
using orbit-based models and using data that extend to 
similar radii, though they do not attempt to characterize 
the shape of the halo. Thus, the comparison of dark 
halo mass results is informative. We find very similar 
numbers. Inside of 39" (which they call Re), they find 
a dark matter fraction of 18%. Inside of 45" (which we 
call Re), we have a dark matter fraction of 13%. There 
are differences in the models as well. First, they use a 
"maximum M/L" model where they force the M/L of 
the stars to have a maximum value. We find the best- 
fit stellar M/L amongst the range modeled. Second, the 
modelling codes are different, with the main difference 
in that they use regularization (which trades the best 
fitted values with smoothness) and we report results for 
the best fit to the data. Third, they use SAURON data 
at large radii and we use our HET data at large radii. 
Their data extend to slightly larger radii (110" compared 
to our limit of 90"), and our data is high signal-to-noise. 
Given all of the these differences, it is impressive that 
we obtain similar results for the dark halo mass. This 
implies that systematic differences are not significant for 
determining the enclosed mass profile. 

5.5. Smoothing 

It is useful to constrain the orbital weighting so that 
the resulting DF is smooth, as a real galaxy's DF may 
be presumed to be. Although we do not usually report 
results when smoothing our models (we argue that al- 
lowing the best fit to the data is the most robust way to 
provide an un-biased result), oth er groups suggest that 
it is important for their model. iRix et all (|1997| ) and 
subsequent studies minimize the variation in the DF, a 
process they term regularization. We employ maximum 
entropy to find the best combinat ion of orbit we i ghts t o 
match the data, as described in I Thomas et al.l J2005) . 
We define a function / = \ 2 — aS where \ 2 is the sum of 
squared residuals to the data, S is the entropy, and a is a 
parameter describing the relative weights of entropy and 
residuals in the fit. In order to minimize / we typically 
start with a large value of a and make it smaller until the 
X 2 no longer varies. To test the effect of smoothing we 
run models such that the it erations stop when a = 0.01, 
a reasonable value based on I Thomas et al.l ()2005j ) . 

Using only our HET data we ran our no-halo and 
NFW-halo models with and without smoothing using a 
coarser grid in parameter space. We find that smooth- 
ing does not alter the results. All of the models have 
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TABLE 5 
Smoothing Model Results 



halo 


smoothing 


x 2 


r s 


c 


P 








kpc 




M Q /pc 3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


none 


no 


940.857 








none 


yes 


1033.73 








NFW 


no 


766.47 


1050 


2.32 


1.14 x 10~ 4 


NFW 


yes 


852.77 


800 


2.70 


1.56 x 10~ 4 



Note. — (1) Dark halo density profile. (2) Smoothing or 
no smoothing. (3) \ 2 of best-fitted model. These \ 2 values 
are lower than those of TableHlbecause the models fit fewer 
data points (HET data only) than the models in Table [4] 
(HET and SAURON). (4) Scale radius of best-fitted model. 
(5) NFW concentration parameter determined from scale 
density. (6) NFW Scale density pcrit^c- 

a lower x 2 using only the HET data than the main re- 
sults of our paper which use both HET and SAURON 
data. The models with smoothing have a larger \ 2 than 
without smoothing (see Table [5]), but the A\ 2 between 
different halo models remains the same. The best-fitted 
NFW dark halo parameters are consistent within the er- 
rors. The internal moments are also consistent with the 
unsmoothed models within the errors. Using an estimate 
by eye, the smoothed model's DF (plotted as I3 versus 
L z in E bins) looks similar to the unsmoothed model's 
DF when smoothed. 

We therefore determine that adding smoothing via 
maximum entropy does not alter the measured halo, in- 
ternal moments, or overall DF shape. We also note that 
these results using only our HET data are consistent with 
those using both HET and SAURON presented through- 
out this paper. 

5.6. Planetary Nebula Data 

We model the NGC 8 21 planetary nebula data of 
iRomanowskv et ail (|2003f) with the best-fitted halos from 
the stellar data. We are not trying to constrain mod- 
els using this data, but rather are interested in what 
orbital properties the planetary nebulae would require 
given the potential derived from the stellar data. In 
doing this we assume that the potential derived from 
stars is correct and that the planetary nebulae are dis- 
tributed in the same w ay as the sta r s. Th is assumption 
may not be realistic, as D ekel et al.l (2005) predicts that 
it is the densities, not the anisotropies, that differ. Fig- 
ure [16] shows the ratio of radial to tangential dispersion 
for the models with no dark halo and best-fitted NFW 
and power- law halos. As expected from the results of 
IRomanowskv et ail (|2003f ). the model with no dark halo 
is roughly isotropic throughout, and tangential at large 
radii. The best-fitted NFW model requires radial or- 
bits throughout and the best-fitted power-law halo re- 
quires extremely radial orbits, with a ra diai / ' a tangential of 
over 3 (corresponding to a /3 of 0.9). This again demon- 
strates the strong mass-anisotropy degeneracy in dynam- 
ical studies. All three models are an excellent fit to the 
data, although there is a preference for a dark halo, but 
it is not statistically significant. 

6. CONCLUSIONS 




— NFW halo 

- -Power law Halo 



1 10 100 

r (arcsoo) 

Fig. 16. — Ratio of radial to tangential dispersion (an average 
of erg and cr^, including streaming motion) along the major axis 
for models with only planetary nebula data (thick lines) and the 
best-fitted no halo, NFW halo, and power-law halo derived from 
the stellar kinematics. All three models are consistent with the PN 
data, in terms of \' 2 , with the dark halo models providing a slightly 
better fit. Thin lines show results for stellar data as comparison. 

We present kinematics of NGC 821 to over 2 effective 
radii using long-slit spectroscopy from the Hobby-Eberly 
Telescope and find that our measured stellar line-of-sight 
velocity distributions are larger than the planet ary neb- 
ulae measurements of lRomanowskv et al.l (|2003f ) at large 
radii. 

Regardless of the density profile used, we are able to 
constrain the enclosed mass of NGC 821 within our kine- 
matic data (~ 2i? e ) as roughly 2 x 1O 11 M0, equally di- 
vided between stars and dark matter. At lR e the ratio 
of dark matter halo mass to total mass is 0.19 for the 
best-fitted NFW halo profile and 0.13 for the best-fitted 
power-law halo profile. 

We find that the best-fitted model of the dark halo in 
NGC 821 has a nearly flat power-law density profile. This 
dark halo gives a better fit than both the NFW halo mod- 
els and models without a dark halo at a greater than 99% 
confidence level. This slope is somewhat unexpected, 
and is strongly inconsistent with halo pr ofiles with inner 
slopes greater tha n one (e.g. isothermal. lHernquistlll99'ol : 
iMoore et al.lll999T ) , and may lend support to halos with a 
flat inner slope (e . g. co red isothermal, logarithmic poten- 
tial, and lBurkertlll995l ). Additionally, one would expect 
that adiabatic co ntraction would create even steeper in- 
ner h alo profiles (jBlumenthal et aD 119861 : iGnedin et al.1 
l200l . which is in conflict with our result. This halo 
result is driven by the data at large radii. 

Our NFW x 2 space shows a degeneracy in radius and 
density as expected. This degeneracy is slightly tilted 
from the expected NFW correlations. Constraining these 
NFW radius and density parameters using a single con- 
centration parameter could lead to biased results. 

In addition to having a significantly poorer fit, the 
models without a dark halo show tangential anisotropy 
at large radii. This may be an indication that a dark 
halo is necessary because the radial component of the 
velocity dispersion may need to be artificially decreased 
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at large radii in order to create a smaller total velocity 
dispersion that can be reproduced by a haloless model. 
The best-fitted dark halo model shows a radial bias in 
the 9 direction at all radii. However we do show that 
the velocities in the <f> direction are greater than the ra- 
dial component. If the planetary nebulae are on radial 
orbits, that would explain why our measured stellar ve- 
locity dispersions ar e larger than the reported planetary 
nebulae dispersions (|Romanowskv et aLll2003|) . We show 
this by modeling the planetary nebula data assuming the 
potential of our best-fitted halo models from the stellar 
kinematics. We find that the planetary nebulae do re- 
quire radially anisotropic orbits to match the best-fitted 
halo potentials. 
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